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Abstract — For fixed compact connected Lie groups H C G, we 
provide a polynomial time algorithm to compute the multiplicity 
of a given irreducible representation of H in the restriction of 
an irreducible representation of G. Our algorithm is based on a 
finite difference formula which makes the multiplicities amenable 
to Barvinok's algorithm for counting integral points in polytopes. 

The Kronecker coefficients of the symmetric group, which 
can be seen to be a special case of such multiplicities, play an 
important role in the geometric complexity theory approach to 
the P vs. NP problem. Whereas their computation is known to be 
#P-hard for Young diagrams with an arbitrary number of rows, 
our algorithm computes them in polynomial time if the number 
of rows is bounded. We complement our work by showing that 
information on the asymptotic growth rates of multiplicities in 
the coordinate rings of orbit closures does not directly lead 
to new complexity-theoretic obstructions beyond what can be 
obtained from the moment polytopes of the orbit closures. Non- 
asymptotic information on the multiplicities, such as provided 
by our algorithm, may therefore be essential in order to find 
obstructions in geometric complexity theory. 

I. Introduction 

The decomposition of Lie group representations into ir- 
reducible sub-representations is a fundamental problem of 
mathematics with a variety of applications to the sciences. 
In atomic and molecular physics (Clebsch-Gordan series), as 
well as in high-energy physics, this problem has been studied 
extensively fl]-[3), perhaps most famously in Ne'eman and 
Gell-Mann's eight-fold way of elementary particles |0]-|6]. In 
pure mathematics, the combinatorial resolution of the problem 
of decomposing tensor products of irreducible representations 
of the unitary group by Knutson and Tao has been a recent 
highlight with a long history of research Q 0. More recently, 
the theories of quantum information ll9l jTTl . computation and 
complexity 0~2), as well as the geometric complexity theory 
approach to the P vs. NP problem |[T3l - [T5l have brought 
the representation theory of Lie groups to the attention of the 
computer science community. 

In this paper, we study the problem of computing multiplic- 
ities of Lie group representations: 

Problem 1.1 (Subgroup Restriction Problem). Let f : H — > G 

be a homomorphism between compact connected Lie groups H 
and G. The subgroup restriction problem for / is to determine 
the multiplicity of the irreducible H -representation Vh,^ 



in the irreducible G-representation Vg,x when given as input 
the highest weights \x and A ( specified as bitstrings containing 
their coordinates with respect to fixed bases of fundamental 
weights, see \§V). 



The name subgroup restriction problem comes from the 
archetypical case where the map / is induced by the inclu- 
sion of a subgroup H C G. IProblem 1. 11 is also known as 
the branching problem. The main result of this paper is a 
polynomial-time algorithm for IProblem 1. 11 

Theorem 1.2. For any homomorphism f : H — > G be- 
tween compact connected Lie groups H and G, there is a 
polynomial -time algorithm for the subgroup restriction prob- 
lem for f. 



Indeed, we describe a concrete algorithm ( [Algorithm V. 1) . 
In particular, for any fixed A and /j the stretching function 
k >->• j7i£* can be evaluated in polynomial time. 

Corollary 1.3. For any homomorphism f : H — > G between 
compact connected Lie groups H and G, positivity of the 
coefficients can be decided in polynomial time. 

Mulmuley conjectures that deciding positivity of the mul- 
tiplicities is possible in polynomial time if the group 



homomorphism / is also part of the input |16|. Corollary 1.3 



can be regarded as supporting evidence that this conjecture 
might in fact be true for general / (note that for specific 
families of homomorphisms, such as those corresponding 
to the Littlewood-Richardson coefficients, positivity can be 
decided in polynomial time [8] El). However, any approach 
to deciding positivity that proceeds by computing the actual 
multiplicities is of course expected to fail, since the latter 
problem is well-known to be ^P-hard IfTSl [19l . 

We establish ITheorem 1. 21 by deriving a novel formula 
for the multiplicities rrC dTheorem IV.2t , which is obtained 
in three steps: First, we restrict from the group G to its 
maximal torus Tq; the corresponding weight multiplicities 
can be computed efficiently by using the classical Kostant 
multiplicity formula 11201 |2D or in fact by evaluating a single 
vector partition function [22-24] ( |§II) . Second, we restrict 
all weights to a maximal torus T# of H. Third, we recover 
the multiplicity of an irreducible iJ-representation by using 



a finite-difference formula (Proposition III. 1 1. By carefully 
combining the first two steps. IProblem 1. 1 1 can be reduced to 
counting integral points in certain rational convex polytopes 
of bounded dimension, which can be done efficiently by using 
Barvinok's algorithm [[2511271 (see also 12811301 ). 

The multiplicity formula itself has intrinsic interest beyond 
its application to algorithmics. One insight that is immediate 
from our result is the piecewise quasi-polynomial nature of 
the multiplicities ( |Corollary IV.3) >. 

Let us now turn to the computation of the Kronecker 
coefficients g\ }ll ^, which arise in the decomposition of tensor 
products of irreducible representations of the symmetric group 
S k GQ: 

where we denote by [A] the irreducible representation of 
Sk labeled by the Young diagram A with k boxes ( |§VI[ ). 
Kronecker coefficients are notoriously difficult to study, and 
finding an appropriately strong combinatorial interpretation is 
one of the outstanding problems of classical representation 
theory. They appear naturally in geometric complexity theory, 
where their efficient computation has been subject to various 
conjectures |[T6l . as well as in quantum information theory in 
the context of the marginal problem and coding theory ifTUl 
Hl][32H3l. 

Using Schur-Weyl duality, the Kronecker coefficients for 
Young diagrams with a bounded number of rows can be equiv- 
alently characterized in terms of a single subgroup restriction 
problem for compact connected Lie groups ( |§ VI) . Therefore, 
by ITheorem 1. 21 they can also be computed efficiently: 

Corollary 1.4. For any fixed d G Z>o, there exists a 
polynomial-time algorithm for computing the Kronecker coef- 
ficient gx.^.v given as input Young diagrams A, /i and v with 
at most d rows. That is, the algorithm runs in (9(poly(log k)) 
where k is the number of boxes of the Young diagrams. 

Corollary 1.5. Positivity of Kronecker coefficients for Young 
diagrams with a bounded number of rows can be decided in 
polynomial time. 

By specializing our technique, we get a clean closed-form 



expression for the Kronecker coefficients (Proposition VI. 1 



which not only nicely illustrates its effectiveness, but also 
implies piecewise quasi-polynomiality for bounded height (a 
feature that has only been noticed in a special case |36|). 
Moreover, it is immediate from our formula that the problem 
of computing Kronecker coefficients with unbounded height 
is in GapP, as first proved in lH9l . 

Similar conclusions can be drawn for the plethysm coef- 
ficients, which can also be formulated in terms of subgroup 
restriction problems [37]. Like the Kronecker coefficients, they 
play a fundamental role in geometric complexity theory lfT5] 
38 1 and quantum information theory ifTTI [39). 

In practice, our algorithms appear to be rather fast as long 
as the rank of the Lie group G is not too large. In the case 



of Kronecker coefficients for Young diagrams with two rows, 
we can easily go up to k — 10 s boxes using commodity 
hardware. In contrast, all other software packages known to the 
authors cannot go beyond only a moderate number of boxes 
(k — 10 2 on the same hardware as used above). Moreover, 
by distributing the computation of weight multiplicities onto 
several processors, we have been able to compute Kronecker 
coefficients for Young diagrams with three rows and k = 10 5 
boxes in a couple of minutesQ We hope that our algorithm will 
provide a useful tool in experimental mathematics, theoretical 
physics, and geometric complexity theory. 

Our final result concerns the asymptotics of multiplicities in 
the general algebro-geometric setup of the geometric complex- 
ity theory approach to proving the VP ^ VNP conjecture, 
an algebraic version of the P ^ NP conjecture. Recall that, in 
a nutshell, this approach amounts to showing that for certain 
pairs of projective subvarieties X and Y one is not contained 
in the other; this would then imply complexity-theoretic lower 
bounds. Both the permanent vs. determinant problem, which 
is equivalent to the VP vs. VNP problem [40|, as well as the 
complexity of matrix multiplication ATI can be formulated in 
this framework IFT3T1T5] l42l . More concretely, let us denote by 
mH,x,k(n) the multiplicity of the dual of an irreducible H- 
representation Vh,^ in the k-th graded part of the coordinate 
ring of X, and similarly for Y (cf. |§VH| for precise definitions). 
Then, 

X C Y => m H ,x,k{tJ) < m H< r, k {p) (11) 

for all ji and k > 0. Therefore, the existence of fi and k such 
that rriH,x,k{li) > mH,Y,k{n) proves that X % Y; such a 
pair (/!, k) is called an obstruction ifPfl . One can relax this 
implication further and instead compare the support of the 
multiplicity functions, 

XCY ^ (m H jcM => m H ,Y,h{n) °) • 

Since computing multiplicities in general coordinate rings is a 
difficult problem, it is natural to instead study their asymptotic 
behavior. Following an idea of Strassen [43 1, it has been 
proposed in l42l to consider the moment polytope, 



A x := 



k=l 



which is a compact convex polytope that represents the asymp- 
totic support of the stretching function. Moment polytopes do 
have a geometric interpretation, which should facilitate their 
computation [44|. Clearly, 



X CY => A x C A 



(1.2) 



However, preliminary results suggest that the right-hand side 
moment polytope Ay might be trivially large in the cases of 
interest |[T5l l42l l45l l46l . and therefore insufficient for finding 
complexity-theoretic obstructions. 

'A preliminary implementation of the algorithm is available upon request 
from the authors. 



It has therefore recently been suggested to study the asymp- 
totic growth of multiplicities (e.g., Il47l §2.2]). The natural 
object is the Duistermaat-Heckman measure, which is defined 
as the weak limit 



DH 



x 



fi/k i 



(1.3) 



where dx <E Z>o is the appropriate exponent such that DHx 
is a non-zero finite measure [4-8 1 . The Duistermaat-Heckman 
measure has a continuous density function fx with respect 
to Lebesgue measure on the moment polytope; it is sup- 
ported on the entire moment polytope (both statements follow 
from the main result of ll48l ). For well-behaved varieties, 
Duistermaat-Heckman measures have a geometric interpreta- 
tion [49-55 1, which makes their computation potentially much 
more tractable ||56H58l (this connection is however less clear 
in the singular cases relevant to geometric complexity theory). 
In this context, our main technical result is the following (see 
[gVHl forthe proof): 

Theorem 1.6. The exponent dx is equal to dim A — Rx, 
where Rx is the number of positive roots of H that are not 
orthogonal to all points of the moment polytope Ax- 

The significance of ITheorem 1. 61 is that the order of 
growth of the "smoothed" multiplicities, as captured by the 
Duistermaat-Heckman measures, does only depend on the 
dimension of the orbit closures and on their moment polytopes. 

Now suppose that we are in the situation that X and Y 
cannot be separated by using moment polytopes, i.e., Ax C 
Ay. For the orbit closures X and Y that one tries to separate 
in geometric complexity theory, one can show that dim A < 
dim!" |Q3] g2). Then, X C Y would imply that d x < dy 
( Corollary VII. 4[ >. But this means that we cannot deduce from 
(I. II) and dl.3t a criterion of the form 



icy 



/x(m)<MaO (Vm), 



since in order to take the weak limit we need to divide 
by different powers of k. Therefore, Duistermaat-Heckman 
measures do not directly give rise to new obstructions, indi- 
cating that a more refined understanding of the behavior of 
multiplicities in coordinate rings might be required. 

II. Preliminaries 

In this paper we will use basic notions of the theory of 
compact Lie groups ll37l I59TI6T1 . Let G be a compact con- 
nected Lie group with Lie algebra q. We fix a maximal torus 
Tg Q G and denote by tc its Lie algebra, the corresponding 
Cartan subalgebra. We write Aq = ker exp | for the integral 
lattice and A G for the weight lattice, which we can consider as 
a subset of t G . The Weyl group Wg acts on t* G by reflections 
through the hyperplanes orthogonal to the roots. Let us choose 
a set of positive roots Rg.+ C A g . This determines a positive 
Weyl chamber t G + , as well as a basis of fundamental weights 



, }, where vq — dimTc is the rank of the Lie 



group, and the Weyl vector P = \ ^Z a 



a. 



The set of 



dominant weights A G , is by definition the intersection of the 
weight lattice and the positive Weyl chamber. 

The fundamental theorem of the representation theory of 
compact connected Lie groups is the fact that the irreducible 
(complex) representations of G can be labeled by their highest 
weight A £ A G + (6TJ; for every element A £ A' G , there ex- 
ists a unique irreducible representation Vg.x with this highest 
weight. Given an arbitrary finite-dimensional (complex) G- 
representation V, we can always decompose it into irreducible 
sub-representations V = © AeA « ma.v (A) Vg,a- We shall 
call the function nicy thus defined the highest weight multi- 
plicity function. 

If we restrict the representation to the maximal torus, we can 
similarly decompose into irreducible representations. Since Tq 
is a compact Abelian group, we can always jointly diagonalize 
its action, and it follows that the irreducible representations 
are one-dimensional. The joint eigenvalues can be encoded 
as a weight f3 6 A G , and we will denote the corresponding 
irreducible representation of T G by C^. The decomposition 
V = (B^gA* m T G ,v{P) C,3 then defines the weight multiplic- 
ity function niT G ,v- We also set [k] = {1, . . . , k}, and write 
/ ~ g for the asymptotic equivalence lim^oo f(k)/g(k) = 1. 

An equivalent way of encoding weight multiplicities is in 
terms of the (formal) character, 



chV 



E 



m T G ,v{P) e p , 



which can be understood as the generating function of rriT G ,v- 
Formally, ch V is an element of the group ring Z[A£.], which 
consists of (finite) linear combinations of basis elements e' 3 
subject to the relation e^e^ = e^ + ^ . The character of an 
irreducible representation Vg,\ is given by the Weyl character 
formula [61, p. 319], 

E w& w G ^We w{X+p) 



chV G .x = 

Observe that we have 
1 



(1 



(II. 1) 



n 



(l-e-°) 



=n< 

a€R G ,. 

= E 



-2o 



(11.2) 



where 4>r g , is the Kostant partition function given by the 
formula 



\ R a,+ \ . ST 



3 



x, = /3}. (11.3) 



That is, <f>R G , counts the number of ways that a weight can 
be written as a sum of positive roots (this number is always 
finite since the positive roots span a proper cone). It follows 
directly from (III. 1 b and (III.2b and that 

/3eA* uiEWg 



chV G: 



A 



In other words, the multiplicity of a weight f3 in an irreducible 
representation Vg,\ is given by the well-known Kostant mul- 
tiplicity formula [20], 



wt g ,v g , a (/3) 



det(w)<f> ROt+ (w(\+p)-p-p). (H.4) 

For any fixed group G, the Kostant partition function can 
be evaluated efficiently by using Barvinok's algorithm [25 1, 
since it amounts to counting points in a convex polytope 
in an ambient space of fixed dimension. Therefore, weight 
multiplicities for fixed groups G can be computed efficiently. 
This idea has been implemented by Cochet lETl to compute 
weight multiplicities for the classical Lie algebras (using the 
method presented in 1 30 1 instead of Barvinok's algorithm). We 
remark that the problem of computing weight multiplicities 
is of course the special case of IProblem 1. 11 where H is the 
maximal torus Tq C G. 

Weight Multiplicities as a Single Partition Function 

If G is semisimple, we can find s,t S Z>o and group 



homomorphisms A: I s 
that 



1} and B : A* 



A 



G 



1 such 



mT G ,Va,M = <t>A{B{l)) (VAeA^^GAy, (H.5) 
where tpA is the vector partition function defined by 



Mv) = #{x g Ko ■■ Ax = y}- 



(11.6) 



Note that this improves over the Kostant multiplicity formula 
( 111.41 . where weight multiplicities are expressed as an alternat- 
ing sum over vector partition functions. In particular, (III. 5b is 
an evidently positive formula. It has been established by Billey, 
Guillemin, and Rassart for the Lie algebra su(d) [22 1, and was 
later extended to the general case by Bliem [23 1 by considering 
Littelmann patterns [62| instead of Gelfand-Tsetlin patterns 



The assumption of semisimplicity for (III.51 > is not a restric- 
tion. Indeed, if G is a general compact connected Lie group 
then its Lie algebra can always decomposed as 



= [0,9] ©3, 



(II.7) 



where the commutator [q, q] is the Lie algebra of a compact 
connected semisimple Lie group G ss , and where 3 the Lie 
algebra of the center Z(G) of G EU Corollary 4.25]. Let 
us choose a maximal torus Tg ss of G ss that is contained in 
Tq. Consider now an irreducible G-representation Vg,\ with 
highest weight A. By Schur's lemma, each element in Z(G) 
acts by a scalar. Therefore, all weights j3 that appear in the 
weight-space decomposition have the same restriction to 3. It 
follows that 



m T G ,V G , A (/3) = 



I otherwise, 



(II.8) 



where we write /i ss and \i z for the restriction of a weight fi 
to the Cartan subalgebra of [q, q] and to 3, respectively. These 
multiplicities can therefore be evaluated by using (III. 5) . 



III. The Finite Difference Formula 

Let V be an arbitrary finite-dimensional representation 
of the compact, connected Lie group G. Clearly, we can 
compute the weight multiplicity function rriT G .v from the 
highest weight multiplicity function may by using any of the 
classical formulas ( III. Il l and (III.4I I. or by evaluating the vector 
partition function ( III. 5b described in |§II| By "inverting" the 
Weyl character formula, the converse can also be achieved: 

Proposition III.l. The highest weight and weight multiplicity 
function of a finite- dimensional G-representation V are related 
by 



m G ,v 




-D a rnr, 



G,V 



where (D a m)(X) = m(A + a) — m(A) is the finite- difference 
operator in direction a. Note that any two of the operators D a 
commute, so that their product is independent of the order of 
multiplication. 

Proof: By linearity, it suffices to establish the lemma for 
a single irreducible representation V = Vg,\ of highest weight 
A. The Weyl character formula ( III. It can be rewritten in the 
form 

w{\+p)-p 



Y det(w)e w ^ +p >- p . (III.l) 



n i i -^ a ) chv cx 

a>0 w£W G 

If we identify elements in Z[AJj] with functions on the weight 
lattice, applying finite-difference operators D a corresponds to 
multiplication by (e~ Q — 1). Therefore, the left-hand side of 

(PCD is identified with (ria e _R G , + -^oWr Q ,Vb, v 

Now consider the right-hand side of (IIII.U . Since X + p is a 
strictly dominant weight, it is sent by any Weyl group element 
w 7^ 1 to the interior of another Weyl chamber. That is, there 
exists a positive root a € Rg,+ suc h that (a, w(X + p)) < 0. 
In particular, w(X + p) — p is never dominant unless w = 1. 
It follows that the restriction of ^Ilae.R G + — ^a) Tnr Gt Va \ 
to Ag , is equal to the indicator function of {A}, i.e., equal 
to the highest weight multiplicity function of Vg,\- ■ 
The idea of using dll.ll ) for determining multiplicities of 
irreducible representations goes back at least to Steinberg 
[64 1, who proved a formula for the multiplicity c v x of an 
irreducible representation Vg,v in the tensor product Vq,x ® 
V(3 )M . These multiplicities c\ are called the Littlewood- 
Richardson coefficients for G. Steinberg's formula involves 
an alternating sum over the Kostant partition function (III. 4b ; 
it can be evaluated efficiently as described by Cochet ETI . 
De Loera and McAllister give another method for computing 
Littlewood-Richardson coefficients ll65l . which applies Barvi- 
nok's algorithm to results by Berenstein and Zelevinsky [66|. 
Since the tensor products of irreducible G-representations are 
just the irreducible representations of G x G, the problem 
of computing Littlewood-Richardson coefficients is again a 
special case of IProblem 1.11 The following consequence of the 



proof of Proposition III.l will be convenient in the sequel: 



Corollary III.2. Write U ae R G , + i 1 - e ~") = S 76 r G ^ 
with Tq C A* g finite and all c 7 =^ 0. Then, 

m G y(\)= ^2 c i m T G ,v{^ + l)- 



In particular, it is evident from |Corollary HL2] that, for any 
fixed group G, the multiplicity of an irreducible representation 
in some representation V can be computed efficiently from 
the weight multiplicities of V by computing a finite linear 
combination. 

IV. Multiplicities for the Subgroup Restriction 
Problem 

Every G-representation V can be considered as ("restricts 
to") a representation of H by setting 



h-v:=f(h)-v {VheH), 



(IV. 1) 



and the subgroup restriction problem for /, as defined in 
Problem I. II amounts to determining the multiplicity of 
a given irreducible representation of H in the restriction of 
a given irreducible representation of G. In this section we 
will derive a formula for these multiplicities (ITheorem IV. 2b . 
which will be the main ingredient of the algorithm presented 
in |§V| below. It will also follow from this formula that the 
are given by a piecewise quasi-polynomial functiorQ in A and 
p ( [Corollary IV3) . 

Let us choose the maximal torus Tjj C H in such a way that 
/{Th) C Tq, and denote the corresponding Cartan subalgebra 
by tff. Of course, this implies that the induced Lie algebra 
homomorphism Lie(/) sends the Cartan subalgebra of H in 
the one of G. Since / is a group homomorphism, Lic(/) 
restricts to a homomorphism between the integral lattices, 
F: A H — > A G , X i-> Lie(/)X. The dual map between the 
weight lattices is given by 



F* : A, 



G 



A 



Hi 



/3h-> (5oF= /3oLie(/)l 



(IV.2) 



The following is well-known and easily follows from the 
definitions: 

Lemma IV. 1. Let V be a representation of G and v £ V a 
weight vector of weight j3 £ A* G . If we restrict the action to H 
via ( II V. Il l then v is a weight vector of weight F*(j3) 6 A H . 

Let us also fix systems of positive roots Rh,+ for H. This 
in turn determines the set of dominant weights A* H + as well 
as a basis of fundamental weights (wj*) as described in |§II| 
Let us also set ru = dimTjf. 

Our strategy for solving the subgroup restriction problem 
for / then is the following: Given an irreducible represen- 
tation Vg,\ of G, we can determine its weight multiplic- 
ities with respect to the maximal torus Tq by using any 
of the formulas presented in |§II| We then obtain weight 
multiplicities for Th by restricting according to ILemma IV. 1 1 

2 In the context of this paper, a quasi-polynomial function is a polynomial 
function with periodic coefficients; see p.|6]for the precise definition. It should 
not to be confused with the notion of quasi-polynomial time complexity. 



Finally, we reconstruct the multiplicity of an irreducible 
representation Vh,h by using the finite-difference formula 
( Proposition III. HCorollary III.2 >. If this procedure was trans- 
lated directly into an algorithm, the runtime would be poly- 
nomial in the coefficients of A (with respect to the basis 
of fundamental weights), i.e., exponential in their bitlength, 
since the number of weights is of the order of the dimen- 
sion of the irreducible representation Vg,x, which according 
to the Weyl dimension formula is given by the polynomial 

n 



aeRc, 



(a, A + p) I (a, p) (cf. the formula by Straumann 
[67 1). We will now show that it is possible to combine the 
weight multiplicity formula (III.5b with the restriction map F* 
in a way that will later give rise to an algorithm that runs in 
polynomial time in the bitlength of the input: 

Theorem IV.2. Let f : H — > G be a homomorphism of com- 
pact connected Lie groups. Then we can find s, s' ,u £ Z>o 
and group homomorphisms A: 1 



s+s 



— > 



and B: A G 



A* H 1 U with the following property: For every irreducible 
representation Vg,\ of G and Vh,h of H, the multiplicity 
of the latter in the former is given by 



c 7 #{x G Z| © I s ' : Ax = B 



p + 1 



where the (finite) set T u and the coefficients (c 7 ) are defined 

b y T\a£R Ht+ C 1 ~ e_Q ) = E 7 er H c 7 e ~ 7 and c 7 °- In fact, 
we can choose s = 0(r G ), s 1 < tq and u = 0(r G ) 



TH- 



Proof: By definition and Corollary III.2| we have 



1 £V H c 7 m T H ,Va, 



(p + 7). In view of 



ILemma IV. II the multiplicity of a Ty-weight 5 6 A^ in the 
irreducible G-representation Vg,\ is given by 



m T H y G A 5 ) 



F*(P)=5 



As in (III. 71 ). let us now decompose the Lie-algebra q — [g, g] 
3. Denote the Lie group corresponding to [g,g] by G ss and 
choose a maximal torus Tg ss which is contained in T. Using 

dira. 



m T G ,V G , x {P) = 



F*(/3)=S 



C BS /3 ss +C a A 



,(A»). 



where we have decomposed F* as a sum of two homomor- 
phisms G ss : Ag aa -> A* H and C z : A* z , G s -> A^. 

Let us now choose group homomorphisms A : Z s —> Z* and 
B = B 1 ®B 2 : A* Gss © A^ s Z* such that $0$ holds for 
the weight multiplicities for G ss . For this, s and t can be taken 



of order 0{r G ) J23] Proposition 19]. Then, 



In particular, this implies that the stretching function k 



'"/<: ,Va .\ (/3ss) 



C B3 /3 sa +C z \z=S 



C HS /3 HS +C,A,=5 



(fc)> 



#{(*,A.):(^ a )(/S.) = (_S^r +a )} 
#{(^As):(r^)(AJ = ( B o 1 -°c/x)(v)}. 




(IV.3) 

After choosing a basis of the lattice A G we arrive at the 
asserted formula (with s' = dbxiTa M and u = t + th)- ■ 

We stress that the proof of ITheorem IV. 2 1 is constructive: 
The maps A and B, whose existence is asserted by the 
theorem, are defined in (II V. 3) in terms of A and B, whose 
construction is described explicitly in [22 Proof of Theorem 
2.1] (for the case of & = su(d)) and in J23 §4] (for the 
general case). See |§VI| for an illustration in the context of the 
Rronecker coefficients. 

If one uses the Kostant multiplicity formula dll-4b instead of 
ill.Si in the proof of lTheorem IV.2l then one arrives at a similar 
formula for the multiplicities involving an additional 
alternating sum over the Weyl group of G. After completion 
of this work, we have learned of |j49l Lemma 3.1] which is 
derived in this spirit. 

Piecewise Quasi-Polynomiality 

Let us use the fundamental weight bases fixed above to 
identify A* G ~ Z r ° and = 1T H . The group homomor- 
phisms A and B correspond to matrices with integer entries, 
which we shall denote by the same symbols. Observe that the 
formula in ITheorem IV. 2 1 in essence amounts to counting the 
number n(y) := # (Aa,b(u) H % s+s ^ of integral points in 
certain rational convex polytopes of the form 



l s+s :xi t ...,x B >0,Ax = By}, 

(IV.4) 

parametrized by y G Z TG+rH . Explicitly, 



^ = c 7 n(A,M + 7)' 



(IV.5) 



It is well-known that n(y) is a piecewise quasi-polynomial 
function in y [68 1. That is, there exists a decomposition of 
jr G +r H mt0 polyhedral chambers such that on each chamber 
C the function n(y) is given by a single quasi-polynomial, 
i.e., there exists a sublattice L C Z TG+TH of finite index 
and polynomials (p z ) with rational coefficients, labeled by the 
finitely many points z G Z ra+TH / L, such that n(y) = P\ y ](y) 
for all y G 17°+™ (cf. [69, §2.2]). We record the following 
immediate consequence: 

Corollary IV.3. For any fixed group homomorphism f : H — >• 
G, the multiplicities are given by a piecewise quasi- 
polynomial function in A and fi. 



,k\ 
'k/i 



is a quasi-polynomial function for large k. This is in fact 
true for all k, as has been observed in [16] (cf. l53l for more 
general quasi-polynomiality results on convex cones, and also 
ll70l for further discussion). 

V. Polynomial-Time Algorithm for the Subgroup 
Restriction Problem 

In this section we will formulate our algorithm for the 
subgroup restriction problem, Problem 1.11 Recall that, by 
(II V.5b . the computation of the multiplicities effectively 
reduces to counting the number of integral points in certain 
rational convex polytopes of the form ( IIV.4) . We shall suppose 
that the highest weights A and /j, which are the input to our 
algorithm, are given in terms of their coordinates with respect 
to the fundamental weight bases fixed in |§IV| Clearly, for 
each of the finitely many 7 G the description of the 
polytope A_4.e(A, fi + 7) (say, in terms of linear inequalities) 
is of polynomial size in the bitlength of the input. It follows 
that Barvinok's algorithm can be used to compute the number 
of integral points in each of these polytopes in polynomial 
time [25 1 (see also 11261 l27l ). This gives rise to the following 
polynomial-time algorithm for Problem 1. 11 thereby establish- 
ing ITheorem 1. 21 

Algorithm V.l. Let f: H — >• G be a homomorphism of com- 
pact connected Lie groups. Given as input two highest weights 
A G A G = 17° and \x G A^ = 1T H , encoded as bitstrings 
containing their coordinates with respect to the fundamental 
weight bases fixed above, the following algorithm computes 
the multiplicity in polynomial time in the bitlength of the 
input: 
m <— 

for all 7 G Th do 

n <— =ff= ^A_4 .b (A, /i + 7) n Z s+S ~\ as computed by 
Barvinok's algorithm (see discussion above) 

m <— m + c^n 
end for 
return m 

Here, Aa.b{v) denotes the rational convex polytope defined 
in (|IV.4| l, and the finite index set V h C A* h as well as the 
coefficients (c 7 ) are defined in the statement of \Theorem IV.2\ 

There are at least two software packages which have 
implemented Barvinok's algorithm, namely LattE [71 1 and 
BARVINOK ||69| |72]. In[|I] we have reported on the perfor- 



mance of our implementation of Algorithm V. 1 for computing 



Kronecker coefficients using the latter package. 

Remark V.2. The existence of a polynomial-time algorithm 
for Problem 1. 1 1 in fact already follows abstractly from 
|Corollary IV.3| since in order to compute we merely have 
to evaluate a fixed piecewise quasi-polynomial function. This 
piecewise quasi-polynomial can be computed algorithmically 
by using a variant of Barvinok's algorithm which is also 
implemented in the BARVINOK package; see ||69l Proposition 
2] and also W\ (5.3.1)]. 



VI. Kronecker Coefficients 

As explained in the introduction, the Kronecker coefficients 
play an important role in geometric complexity theory and 
quantum information theory. In this section, we will describe 
precisely how they can be computed using our methods. 

Let us recall the language of Young diagrams which is 
commonly used in this context [31]. A Young diagram with 
r rows and k boxes is given by an ordered list of integers 
Ai > . . . > A r > with \ = k. It can be visualized as 
an arrangement of k boxes in r rows with Xj boxes in the 
j-th row. We set Xj — for all j > r. We will now consider 
the unitary group U(d), which consists of the unitary d x d- 
matrices. Let us fix a system of positive roots and denote 
the corresponding basis of fundamental weights by (uij). To 
each Young diagram A with at most d rows we associate 
the irreducible representation of U(d) with highest weight 
equal to Y2j=i (^j ~ Every polynomial irreducible 

representation of U(<f) arises in this way. By a slight abuse of 
notation, we identify Young diagrams with the corresponding 
highest weights. More generally, we can associate to every 
integer vector /3 G 1> d the weight X^=i iPi ~ Pj+l) °- 'j> where 
we set fid+i = 0. This defines a bijection between Z, d and 
the weight lattice A^j^ of U(d). In particular, the positive 
roots fixed above correspond to the integer vectors of the form 
(...,0,1,0,. ..,0,-1,0,...). 

The Kronecker coefficient g\ : ^ :V associated with triples of 
Young diagrams A, /i and v with k boxes each and at most a, 
b and c rows, respectively, can then be defined in terms of the 
following subgroup restriction problem of compact, connected 
Lie groups: Let H = U(a) x U(6) x U(c) and G = \J(abc) 
and consider the homomorphism / : H — >• G given by sending 
a triple of unitaries (£7, V, W) to their tensor product U ® 
V W. The Kronecker coefficient g\ jflj v is then given by the 
multiplicity of the irreducible ^-representation Vh,(\,ij,,v) = 
Vu(a),\ ® V\J(b),iJ, ® ^U(c),i/ m the restriction of the symmetric 
power Sym fe (C abc ), which is the irreducible G-representation 
labeled by the Young diagram (fc) consisting of a single row 
with k boxes. That is, 



3a„ 



,(fe) 



(VI. 1) 



This definition in fact does not depend on the concrete values 
chosen for a, b and c, as can be seen by rephrasing it in terms 
of the representation theory of the symmetric group Sk [15, 
§8] (but of course a, b and c have to be chosen at least as large 
as the number of rows of the Young diagrams). Moreover, it 
is evident that the Kronecker coefficients are symmetric in the 
variables A, \i, and v. 

It follows that, for any fixed choice of a, b and c, 



Algorithm V. 1 can be used to compute the Kronecker coef- 
ficient dVI.ll ) given Young diagrams with at most a, b and 
c rows, respectively, in polynomial time in the input size, or 
equivalently in time (9(poly(log k)), where k is the number of 



boxes of the Young diagrams. This establishes Corollary 1.4 



coefficients is known to be #P-hard in general lfT9l ; hence 
we do not expect that there exists a polynomial-time algorithm 
without any assumption on the number of rows of the Young 
diagrams. 

When computing Kronecker coefficients using the above 
method, we are only interested in the representation 
V\j(abc).(k) = Sym fe (C obc ), not in arbitrary irreducible rep- 
resentations of U(abc). By specializing the construction de- 
scribed in ITheorem IV.2I to this one-parameter family of 
representations, we obtain the following result: 

Proposition VI.l. The multiplicity of a weight 8 = 
(S A ,S B ,S C ) S Z a © Z b © Z c = A|j (we use the identifi- 
cations fixed at the beginning of \§Vty in the irreducible G- 
representation Sym /c (C af,c ) is equal to the number of integral 
points in the rational convex polytope 



A(k,5) = |(^, m ,„) 



k. 



m.n 



l,m,n 

Ln l.m 



It follows that the Kronecker coefficient for Young diagrams 
A, /i, v with k boxes and at most a, b and c rows, respectively, 
is given by the formula 



3a„ 



E 



c 7 #(A(fc,(A,/^)+7)nZ a6c ) 



where Th and (c 7 ) are defined as in the statement of 
\Corollary lll.2\ 

Proof: It is well-known that the weight spaces for the 
action of U(d) on Sym fe (C d ) are all one-dimensional and that 
the set of weights corresponds to the integer vectors in the 
standard simplex rescaled by k BP . In our case, d = abc, so 
that the weights are just the integral points of the polytope 



vabc 
^>0 



E 



XL 



k 



Moreover, the dual map F* : Ay^ Qfjc j 
defined in (IIV.2I ) is given by 



^lI(a)xU(l))xU(c) aS 



f 1 abc -> 1 a ®z b ®z c 

[( x l,m,n) t— ^ ( X/rn,n x l,m,m x l,rn,n, X/Z,m x l,m 

We conclude that the multiplicity of a weight S — (5 A ,5 B ,S C ) 
for U(a) x U(6) x U(c) is given by the number of integral 
points in the polytope A(k,8) described above. ■ 



Let us again stress that the problem of computing Kronecker 



Just as for our main algorithm, Proposition VI.l gives 
rise to a polynomial-time algorithm for computing Kronecker 
coefficients with a bounded number of rows. This second 
algorithm runs faster than the generic one presented earlier, 
since the ambient space M. abc has a smaller dimension than 
what we would get from the construction described in the 
proof of ITheorem IV.2I We remark that the time complexity 
for unbounded a, b and c can be deduced from 11271 . 



VII. ASYMPTOTICS 

In this section we will prove our result on the generic order 
of growth of multiplicities in the coordinate ring of a projective 
variety ( ITheorem I.6I >. 

We will work in the following general setup: Let V be a 
finite-dimensional rational representation of H, and suppose 
that A is an iJ-stable closed subvariety of the associated 
projective space P(V). The homogeneous coordinate ring 
C[A] is graded, and we can decompose each part into its 
irreducible components, 

oo oo 

c[x\ = 0qi] t = ®®m ttl M Vh.^ (vn.i) 



for the representations that occur in C[A]. The coefficient 

p (v) = UaeRH.+ .alAx (<*>/*>/(<*> P) is a polynomial func- 
tion in fj,. Since Ax is compact, we can therefore find a 
constant C > such that 



dmxV klM < Ck Rx 
It follows that 



(V^eA x n-Ay. 



m H ,x,k(kn) dimVfc^ 



ck Rx+d x / dDH 



fc=0 



fc=0 M 



where, following the usual conventions, we have decomposed 
with respect to the dual representations . The stretching 
function is then by definition k h4 rriH,x,k{kfJ>). We stress that 
in contrast to [16], where it was assumed that X has at most 
rational singularities, we do not even require that X is a normal 
variety ll73l . This is highly relevant for geometric complexity 
theory, since it was recently shown in [74] and [42] that the 
studied varieties (the orbit closures of the determinant and 
permanent on the one hand, and of the matrix multiplication 
tensor and the unit tensor on the other hand) are in fact never 
normal except in trivial situations. 

Remark VII. 1. The subgroup restriction problem for a rational 
group homomorphism / : H — > G can be realized in the above 
setup: Indeed, for any highest weight A £ A G , consider 
X = Og,\, the coadjoint orbit through A, with the induced 
action of H. This variety can be canonically embedded into 
projective space as the orbit of the highest weight vector in 
¥(Vg,x), and it is a consequence of the Borel-Weil theorem 
that C[O a ,x] = (BT=o V G.k\- fi y comparing with dVTTTT l it 

follows that TTlH Or- 



< Ck Rx ^TO# : x,fc(M 

so that dim X < Rx + dx ■ 

On the other hand, since DHx is Lebesgue-absolutely con- 
tinuous, the boundary of the moment polytope does not carry 
any measure. We can therefore find a compact set K contained 
in the (relative) interior of the moment polytope which has 
positive measure with respect to DHx. Note that P(pi) is 
positive for all /z contained in the interior of the moment 
polytope (indeed, for all positive roots a with a JL Ax there 
exists v e Ax such that (a, v) > 0; since we can always write 
fj, as a proper convex combination of v and some other point 
v 1 G Ax, it follows that (a, /i) > 0). This implies that on the 
compact set K we can bound P(/x) from below by a positive 
constant. Thus there exists a constant D > (depending on 
K) such that 

1 



(V/xeXn-A^ + ). 



dim V klJi > D k Rx 
Consequently, 

^2mH,x,k(kfJ.) dimVfc^ > ^m HiX M( k ^) dim ^V 



/xeifniAjj 



,fcO) = In particular, the above > Dk Rx V] m ff) x,* ~Dk Rx+dx [ dDH x 

:ching function, k i-> m H ,o G x ,k(ku), „ c .„i.. Jk 



definition of the stretching function, k (->• rriH,o G 
coincides with our previous usage, k m kx 



k // * 



Proof of \Theorem I.6\ By the Hilbert-Serre theorem, the 
function k dimC[A]fe is a polynomial of degree dim A 
for large k ll73l Theorem 1.7.5]. Hence there exists a constant 
A > such that 

Ak dimX ~dimC[A] fe 
= y^mff.x.fcQ-Q dimV^ = ^ rn H ,x,k{kjj) dimV fcM , 

where for the last equality we have used the definition of the 
moment polytope Ax- By the Weyl dimension formula, we 
have 



We conclude that also dim A > Rx 
equality. 



dx, hence we have 



dimVfc^ = Yl 

1 , ^ 

R 

\ a/Ax / 



(a, kfJL + p) 



k Rx +0{k Rx ~ 1 ) 



Let us now elaborate on the argument presented at the end of 
the introduction, where we showed that Duistermaat-Heckman 
measures do not directly give rise to new complexity-theoretic 
obstructions. For this, we consider a pair of projective sub- 
varieties A and Y with dim A < dimF, as is the case 
for the orbit closures of relevance to GCT Let us assume 
that Ax C Ay, so that the moment polytopes alone do not 
already give rise to an obstruction. Clearly, this implies that 
Rx < Ry- 

Lemma VII.2. Let Ax C Ay and Rx < Ry- Then, 
dim Ax < dim Ay. 

Proof: Note that we have 

dim Ax = dim aff Ax < dim aff Ay = dim Ay , 

with equality if and only if the two affine hulls aff Ax Q 
aff Ay are equal. 



Now by assumption there exists a positive root a <E Rh,+ 
that is orthogonal to all points in Ax (i.e., for all p G Ax, 
a _L p), but not to all points in Ay. It follows that a is 
also orthogonal to all points in the affine hull of Ax, but 
not to all points in the affine hull of Ay. Therefore, we have 
affAxCaffAy. ■ 

Lemma VII.3. Let dim Ax < dim Ay. Then, X CY implies 
dx < dy. 

Proof: If X C Y then it is immediate from (II.lt and 
( 11.31 ) that dx < dy- Let us suppose for a moment that in fact 
dx = dy. Then it follows from ( II. lb that 

f dBRx(ii)g^) < f dDHy (//)<?(//) 

J A x J A Y 

for any test function g. In particular, this inequality would 
hold for g the indicator function of Ax- But this is clearly 
impossible, since DHy is absolutely continuous with respect 
to Lebesgue measure on Ay, for which Ax is a set of measure 
zero. ■ 

Corollary VII.4. Let dim X < dim Y. Then, X CY implies 

dx < dy. 

Proof: Clearly, X C Y implies that Ax C Ay and 
Rx < Ry- If Rx = Ry then the assertion follows directly 
from lTheorem I.6l since 

dx = dim X - R x < dim Y - Ry = dy. 

Otherwise, if Rx < Ry, it follows from combining 
ILemma VTI21 and ILemma VII31 ■ 
As described in the introduction, the upshot of the above is 
that we cannot directly deduce from (II. Il l a new criterion for 
obstructions based on the Duistermaat-Heckman measure that 
goes beyond what is provided by the moment polytope. 
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